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Detection of the fluctuations in 21 cm line emission from neutral hydrogen during the 
Epoch of Reionization in thousand hour integrations poses stringent requirements on 
calibration and image quality, both of which necessitate accurate primary beam models. 

The Murchison Widefield Array (MWA) uses phased array antenna elements which 
maximize collecting area at the cost of complexity. To quantify their performance, we 
have developed a novel beam measurement system using the 137 MHz ORBCOMM 
satellite constellation and a reference dipole antenna. Using power ratio measurements, we 
measure the in situ beampattern of the MWA antenna tile relative to that of the reference 
antenna, canceling the variation of satellite flux or polarization with time. We employ 
angular averaging to mitigate multipath effects (ground scattering), and assess 
environmental systematics with a null experiment in which the MWA tile is replaced with 
a second reference dipole. We achieve beam measurements over 30 dB dynamic range in 
beam sensitivity over a large field of view (65% of the visible sky), far wider and deeper 
than drift scans through astronomical sources allow. We verify an analytic model of the 
MWA tile at this frequency within a few percent statistical scatter within the full width at 
half maximum. Towards the edges of the main lobe and in the sidelobes, we measure tens 
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of percent systematic deviations. We compare these errors with those expected from 
known beamforming errors. 


1. Introduction 

The prospects of studying the formation of the 
first structures in the universe at z ^ 6 and ear¬ 
lier with 21 cm hydrogen emission have driven in¬ 
vestment in a new generation of low frequency ra¬ 
dio astronomy instruments (see Furlanetto et al. 
[2006]; Morales and Wyithe [2010]; Pritchard and 
Loeb [2012]; Loeb and Furlanetto [2013]; Zaroubi 
[2013] for reviews). The extreme surface brightness 
sensitivity required to detect the 21 cm signal in the 
presence of galactic, extragalactic, and thermal noise 
backgrounds has pushed this generation of experi¬ 
ments into an untested regime. Uncertain primary 
beams and source catalogs, in addition to wide fields 
of view, complicate analysis and demand new meth¬ 
ods of calibration, imaging, and primary beam char¬ 
acterization. 

A number of experiments now operating, such as 
the Murchison Widefield Array (MWA) [Lonsdale 
et al, 2009; Tingay et al, 2013; Bowman et al, 

2013] , the Donald C. Backer Precision Array for 
Probing the Epoch of Reionization (PAPER) [Par¬ 
sons et al, 2010, 2014], and the LOw Frequency Ar¬ 
ray (LOFAR) [van Haarlem et al, 2013], as well 
as demonstrator instruments like the MIT Epoch 
of Reionization Array (MITEoR) [Zheng et al, 

2014] , next-generation experiments like the Hydro¬ 
gen Epoch of Reionization Array (HERA; Pober 
et al. [2014]; http://reionization.org), and future in¬ 
struments like the Square Kilometer Array [Hall, 
2005] have opted for large arrays (hundreds of ele¬ 
ments) of non-pointing or only coarsely-pointing an¬ 
tennas, attempting to balance collecting area and 
cost considerations. For all these experiments, in situ 
high-fidelity primary beam characterization remains 
a major challenge given the high dynamic range [e.g. 
Parsons et al, 2012; Beardsley et al, 2013; Thya- 
garajan et al, 2013; Pober et al, 2014] thought nec¬ 
essary to reveal the cosmological 21 cm signal. 

Recovering the 21 cm signal in the presence of 
strong foregrounds is made easier by taking advan¬ 
tage of an effective containment of smooth spectrum 
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foregrounds in a Fourier space region knows as the 
“wedge”, despite the frequency-dependent response 
of the interferometer [Datta et al, 2010; Dillon et al, 
2014a; Pober et al, 2013; Morales et al, 2012; Vedan- 
tham et al, 2012; Thyagarajan et al, 2013; Trott 
et al, 2012; Liu et al, 2014a, bj. However, frequency 
dependent systematics due to insufficiently accurate 
primary beam modeling for calibration or primary 
beam correction may cause foreground leakage out 
of this compact region. This would shrink the region 
within which a cosmological power spectrum mea¬ 
surement can be made and thus lower the significance 
of a detection [Pober et al, 2014]. In fact, Thyagara¬ 
jan et al [2015] show that most pernicious for such 
measurements is sky emission from large zenith an¬ 
gles, even near the horizon, just where beampatterns 
are most difficult to model. Ultimately, the full po¬ 
larization response may be needed to best model and 
subtract polarized sources that can leak a sinusoidal 
frequency signal into Stokes I due to their Faraday ro¬ 
tation [Jelic et al, 2010; Moore et al, 2013]. Recent 
measurements [Bernardi et al, 2013; Moore et al, 
2013; Moore et al, 2015; Asad et al, 2015] indicate 
that most of the high rotation measure sources are, 
however, largely depolarized at low frequency. 

Though brought to the forefront again by 21 cm 
science, primary beam measurements have a long his¬ 
tory in radio astronomy and electrical engineering. 
Radio astronomers typically rely upon celestial radio 
sources with known flux densities to measure beam- 
patterns as the sources trace out cuts through the 
beam [e.g., Thyagarajan et al, 2011]. If the beam is 
narrow enough for the sky to appear as a single point 
source, knowledge of its flux density is not needed 
to measure relative beam sensitivity along its track 
[e.g. Colegate et al, 2014], though combining tracks 
from different sources, or using fields with multiple 
sources requires accurate knowledge of their relative 
fluxes. Indeed, the wide fields of view of dipole ele¬ 
ments and uncertainties in low frequency source cat¬ 
alogs make this analysis difficult and entangled with 
calibration [Jacobs et al, 2013]. Further, the lack 
of axial symmetry in non-dish antennas around the 
antenna pointing direction makes a complete beam- 
pattern impossible to measure from just a handful of 
cuts. Relying only on the weaker assumption of 180° 
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rotation symmetry, Pober et al. [2012] present an in¬ 
terferometric beam measurement technique making 
use of celestial sources with unknown flux densities, 
assuming the data are already calibrated. 

In this paper we pursue an alternate approach 
based on probe signals from satellites. Satellite- 
based antenna beampattern measurements have 
many clear advantages over astronomical sources. 
Satellites are substantially brighter, and thus dom¬ 
inate in otherwise crowded fields and probe deep 
into beam sidelobes. They also make many cuts 
through the beam over the course of many orbits, due 
to precession. Brueckmann [1963] exploited these 
advantages in early satellite measurements of an¬ 
tenna beampatterns. Lingering issues such as vary¬ 
ing antenna pointing and plane of polarization (due 
to Faraday Rotation or simple projection effects) 
and time-varying satellite transmitter power may be 
solved through use of a simple, well-understood refer¬ 
ence antenna and power ratio measurements, which 
measure the relative beampattern of the Antenna- 
Under-Test (AUT) [eg. Fukao et al., 1985; Law et al., 
1997; Hurtado et al, 2001; Butler, 1998]. This ap¬ 
proach has been used to great effect in holographic 
antenna measurements [Rochblatt and Seidel, 1992; 
Harp et al., 2011; Lasenby, 1985; Godwin and et al., 
1986; Deguchi and et al., 1993], with at most a hand¬ 
ful of satellites. 

As a first step towards antenna beam measure¬ 
ments exploiting all these advantages for 21 cm cos¬ 
mology and for the MWA, we develop a prototype 
of a novel beam measurement system using probe 
signals from the 137 MHz ORBCOMM satellite con¬ 
stellation and a reference dipole antenna. The pre¬ 
cession of these low Earth orbit satellites and their 
sheer number (~ 30) yield 65% coverage of the vis¬ 
ible sky (limited by satellite coverage at our Green 
Bank site) at 2° resolution in a single day. Early 
tests and demonstrations of our initial concept were 
presented by Ries [2007]; Bradley and Ries [2008]; 
Czekala and Bradley [2010a, b]; Neben et al. [2014], 
and recently, Zheng et al. [2014] have demonstrated 
a simple implementation of this concept in a working 
interferometer. We use an ORBCOMM interface box 
to determine satellite transmission frequencies on the 
fly, allowing us to take advantage of even more of the 
ORBCOMM signals in an automated manner. 

In this work, we present the full working power 
pattern mapping system, as well as an error analy¬ 
sis of environmental systematics such as multipath 


(ground scattering). We expect the lessons learned 
about these systematics to inform future beam mea¬ 
surements and array calibration methods employing 
probe signals carried by remote-controlled drones or 
satellites. 

We use our working beam measurement system to 
make the first precision measurements of an MWA 
tile in a deployment-style environment over a large 
fraction of sky. Each MWA tile is a 4 x 4 grid of 
bowtie dipoles, optimized to have a broad frequency 
response over 80-300 MHz, whose signals are com¬ 
bined in a delay-line beamformer. This design results 
in a large collecting area per tile and a beam narrow 
enough to steer away from the bright galactic disk 
and terrestrial RFI, but adds model complexity and 
uncertainty near the edges of the main beam and 
in the sidelobes (grating lobes). Indeed the MWA 
tile design poses simulation challenges due both to 
potential dipole cross coupling effects as well as the 
large number of degrees of freedom which must by 
simulated (hundreds of frequencies, dozens of point¬ 
ings, and fine angular resolution). Note, though, 
that current processing of MWA data utilizes sim¬ 
ulated beampatterns, and experiments such as that 
presented in this work aim more to assess their va¬ 
lidity than to replace them with measurements. 

Early beampattern measurements of deployed pro¬ 
totype MWA tiles using source drift scans [Bowman 
et al., 2007] and anechoic chambers [Williams, 2012] 
revealed rough agreement with models. Yet ~ 1 
dB deviations were observed throughout the main 
lobe and 3 — 5 dB deviations in sidelobes, highlight¬ 
ing the need both for better modeling and better 
control of measurement systematics. Later, Sutinjo 
et al. [2015] found that upwards of 200 MHz, inter¬ 
actions between antennas necessitate more complex 
modeling than simple multiplication of a Hertzian 
dipole beampattern by the array factor. We are in¬ 
terested here in characterizing the deviations from 
ideality at lower frequencies where the simple model 
is more likely to hold. At some level, the beams 
are expected to be corrupted by beamforming am¬ 
plitude and phase errors, finite ground screen effects, 
as well as any dipole-dipole interactions. We com¬ 
pare our measurements with an empirical budget of 
beamforming errors due to dipole phase and gain 
mismatching which we present in a separate paper 
(Neben et al., in prep). 

In Section 2 we discuss our beam measurement 
system including the ORBCOMM satellites, our ref¬ 
erence antenna, and our data acquisition system. We 
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discuss our data analysis pipeline in Section 3 and 
also present results of a null experiment in which 
a second reference antenna is used as the Antenna- 
Under-Test. We present beam measurements of our 
MWA tile and compare with models in Section 4, and 
conclude with discussion and conclusions in Section 
5. 

2. Measurement System 

2.1. Overview 

Figure 1 shows a schematic diagram of our antenna 
setup and measurement chain. A reference antenna 
(see Section 2.3) and an Antenna-Under-Test (AUT) 
both pick up transmitted signals from a passing OR¬ 
BCOMM satellite, which are mixed down, sampled, 
and fast fourier transformed by our data acquisition 
system (see Section 2.4), and finally the power as a 
function of frequency from both antennas is saved to 
disk at every time step. Note that in practice, both 
our antennas are dual-polarization, necessitating two 
additional measurement chains like these. Addition¬ 
ally, an ORBCOMM Interface Box (typically sup¬ 
plied to commercial users of the ORBCOMM con¬ 
stellation) interfaces with each passing satellite and 
outputs its identifier (which allows precise prediction 
of the satellite’s location using orbital data) and fre¬ 
quency channel (which specifies which ~20 kHz wide 
ORBCOMM band between 137-138 MHz the satel¬ 
lite is transmitting on). 

This work was performed at the National Radio 
Astronomy Observatory site in Green Bank, West 
Virginia, located in the US National Radio Quiet 
Zone. Though the ORBCOMM satellites domi¬ 
nate the radio sky in the 137-138 MHz band, wide 
band interference from terrestrial radio transmitters 
causes saturation problems in other geographic lo¬ 
cations. Our AUT and reference antenna are posi¬ 
tioned 50m apart on a North-South baseline located 
at (38.429348°N, -79.845737°E), and aligned to an 
accuracy of about a degree, as confirmed by Google 
Earth (http://www.google.com/earth) imagery. The 
Green Bank site is not perfect, however, and the sur¬ 
rounding hills and radio telescopes raise concerns of 
shadowing and multipath effects. We assess these 
with a null experiment in which the beampattern of 
a known reference-style dipole is used as the AUT 
(see Section 3.3). 

In this work we measure the angular response of 
each instrumental polarization (NS and EW) to un¬ 


polarized radiation, leaving for future work measure¬ 
ment of the full polarized beampatterns. While it 
is not obvious that fully polarized satellite probe sig¬ 
nals suffice to measure unpolarized beampatterns, we 
show in Appendix A that this is in fact possible if 
both our reference antenna and our AUT have the 
same polarization response. 

2.2. ORBCOMM Satellite Constellation 

ORBCOMM Inc. operates a constellation of ~ 30 
communications satellites in low Earth orbit (alti¬ 
tude ~ 800km) designed for users requiring low baud 
rate communication with remote sites. The satellites 
provide excellent Earth-coverage and near continu¬ 
ous transmission, predominantly occupying orbital 
planes with inclinations between ±45° and 10 nar¬ 
row (~ 20 kHz wide) subbands in the 137-138 MHz 
band. An advantage of these satellites over higher al¬ 
titude satellites (such as the GPS constellation) for 
beampattern measurements is that good sky cover¬ 
age is achieved far more quickly due to the shorter 
orbital period and rapid orbital precession resulting 
from the lower altitude. In particular, sky coverage 
at our Green Bank site is limited only by absence of 
satellites with inclinations greater than 45°. The in¬ 
formation content of the transmitted satellite signals 
is irrelevant for our purpose and is lost in the RMS 
power measurements of our data acquisition system 
(see Section 2.4). 

Each satellite’s frequency is relatively stable over 
days, but shifts periodically to avoid interference 
within the constellation. There are typically several 
ORBCOMM satellites above the horizon at any given 
time, and while we can easily compute their posi¬ 
tions using published orbital elements, we must de¬ 
termine which frequencies correspond to which satel¬ 
lites. Zheng et al. [2014] use interferometric phases 
to identify and exclude times when more than one 
satellite is present. We are able to take advantage of 
all satellite passes using an ORBCOMM User Inter¬ 
face Box (typically provided by ORBCOMM Inc. to 
users) connected to a separate antenna, whose debug 
port logs the satellite number and frequency band 
occupied by each passing satellite. 

During data collection, we record the satellite po¬ 
sitions using the Linux program predict ^ which nu¬ 
merically integrates the orbits using orbital elements 
(TLE) data published by Celestrak^ daily. We run 
two copies of predict in live multi-satellite track¬ 
ing mode on our data acquisition (DAQ) computer, 
and query them for the angular positions of all satel- 
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lites currently above the horizon whenever a satellite 
power measurement is made (typically every 200ms). 
We save this information with the recorded satellite 
power data. See Section 2.4 for a detailed description 
of our data acquisition system. 

2.3. Reference Antenna 

Our reference antenna is a simple dual-polarization 
dipole mounted above a 2 m x 2 m ground plane, 
elevated 48.3 cm above the soil (see left panel of Fig¬ 
ure 2). The antenna itself is made from 90.4 cm long 
and 0.5 inch diameter copper tubing, and is encapsu¬ 
lated in 2 in diameter Schedule 40 PVC tubing. The 
beampattern of the dipole was derived from an elec¬ 
tromagnetic analysis of the physical structure over a 
finite ground plane using CST Microwave Studio 
The soil is modeled as a cube of lossy dielectric 3.05 
m on a side having a relative permittivity of 13 and 
electrical conductivity of 0.005 S/m. 

2.4. Data Acquisition Hardware 

Our receiver produces a mean power measurement 
in each frequency channel in each of two instrument 
polarizations (NS and EW) for each of our two an¬ 
tennas (the AUT and the reference antenna) every 
~200 ms. The raw signals are first mixed with a 
127 MHz local oscillator down to an intermediate 
frequency of ~10.5 MHz, then simultaneously sam¬ 
pled at 2 MHz, digitally mixing them down to the 
0-1 MHz baseband (i.e. one of the Nyquist zones 
falls in the frequency range 0-1 MHz). A 12-bit 
ADC acquires a burst of 51,200 samples on which 
50 FFTs of length 1024 samples are performed, ef¬ 
fectively covering the 137-138 MHz bandwidth with 
a resolution 2 kHz (cf. ^20 kHz bandwidths of ORB- 
COMMS). From these measured voltage Fourier am¬ 
plitudes Vant,poi,i(/), where i runs from 1 to 50, the 
RMS powers (|14nt,poi,i(/)P)i are saved to disk for 
each antenna, polarization, and frequency channel 
along with a list of all satellites currently above the 
horizon and their angular positions (obtained from 
predict, see Section 2.2). A full complex polar¬ 
ization analysis of the AUT beam would be pos¬ 
sible given measurements of voltage cross powers 
(Rant.poi.*(/)foant'.poi', *(/))*. though this places strin¬ 
gent requirements on instrumental phase stability 
and we thus reserve it for future work. 


3. Data Analysis 

3.1. Processing Satellite Passes 

The first step is to determine each satellite’s trans¬ 
mission frequency during each pass. We have devel¬ 
oped a script to extract this information from the 
captured debug output of the ORBCOMM User In¬ 
terface Box, which periodically “syncs” with passing 
satellites and logs their identifiers and transmission 
frequencies. We assume a time window of 30 min¬ 
utes over which the recorded frequencies are valid, 
and use this mapping of satellites and time windows 
in the next step of the analysis. 

Next, we manipulate the satellite power data. At 
each ~200 ms time step, a beam measurement can 
be made at the position of each satellite above the 
horizon (each of which is transmitting on a differ¬ 
ent, but now known frequency). This is done sep¬ 
arately for EW and NS instrumental polarizations. 
The measured AUT and reference dipole powers in 
each satellite’s frequency band are determined by in¬ 
tegrating the measured RMS band powers over the 
central 15 kHz of the satellite’s signal band. The 
background power level is estimated as the minimum 
of the received power level at start and end of each 
pass (when the satellite is below the horizon), and 
data within 20 dB of that floor are rejected. This 
ensures that the bias on measured beampattern due 
to the sky noise is less than 1%. 

As a heuristic description of our measurements as¬ 
suming the satellite signals are unpolarized, consider 
the received powers by the AUT and the reference 
antenna, Taut and Pref) in one instrumental polar¬ 
ization, and let Href and Paut be their unpolarized 
beam responses at the angular position of the satel¬ 
lite. Let F be the incident flux from the satellite so 
that Pref = Bj-efF and Paut = PautP- Then it is 
clear than the AUT beam response in the satellite 
direction is given by 

Paut = Paut Pref/Pref (1) 

This is essentially our analysis, done separately for 
EW and NS oriented antennas. We show in Ap¬ 
pendix A that the polarization of satellite signals 
does not affect Equation 1, assuming both antennas 
have the same polarization response. 

3.2. Forming Power Patterns 

To form power patterns, we grid the measured 
AUT beam values from many satellite passes into 
equal solid angle cells on the sky using the HEALPix 
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software package [Gorski et al, 2005] to (1) facilitate 
comparison with model power patterns; (2) average 
over short-timescale fluctuations due to noise and 
multipath effects; and (3) facilitate rejection of out¬ 
liers in each sky pixel due to rare saturation issues. 
We choose a cell size of 1.8 deg (nside=32) which re¬ 
sults in ~5 satellite passes and ~ 75 measurements 
in each pixel per day, out of which ~ 5 typically fall 
outside of the central 90% and are rejected as out¬ 
liers. This results in few percent precision on our 
measured beams and sufficient resolution to resolve 
features of interest except within several degrees of 
the the MWA beam nulls, where the beam changes 
on sub-pixel scales. Normalizing the measured beam 
might be done by rescaling it to peak at unity, though 
we opt for a less noisy normalization by fitting for a 
rescale factor to best match the measured beam to 
the normalized analytic model within ^ 10 degrees 
of boresight. 

3.3. Assessing Systematics with a Null Experiment 

We characterize systematics using a null experi¬ 
ment, in which we use a second reference-style an¬ 
tenna as the AUT. The beampatterns of the two ref¬ 
erence antennas will deviate from each other due to 
environmental effects (e.g., multipath effects or shad¬ 
owing) and instrumental non-idealities (e.g., align¬ 
ment errors or imperfections in the ground screen, 
soil, or dipole itself). Our null experiment will effec¬ 
tively measure the ratio of these two antenna beam- 
patterns, and thus, the level at which they deviate 
from each other. We interpret this measure as a 
rough proxy for deviation of each antenna away from 
the ideal electromagnetic model. 

Figure 3 shows a satellite pass from this setup in 
detail. Over the course of 15 minutes, the satellite 
rises out of the background, passes through the vis¬ 
ible sky, and falls below again (shown in the top 
panel). The units are dB relative to the background 
level (shown below to be predominately diffuse galac¬ 
tic emission), and we mark with vertical lines the re¬ 
gion within which the received power is more than 
20 dB stronger than the background. If the AUT 
beampattern is different from that of the reference 
dipole, beam measurements outside this region will 
suffer systematics at a level of a few percent and 
larger due to the diffuse galactic emission received in 
addition to that from the satellite. We opt to simply 
avoid this region in lieu of subtracting a background 
model. The ratio of the two reference antenna pow¬ 


ers (shown in the middle panel) is mostly consistent 
with unity up to 0.5 — 1 dB statistical scatter, sys¬ 
tematic biases of comparable magnitudes are appar¬ 
ent at large zenith angles. As this is just one satellite 
pass, it is difficult to draw general conclusions about 
which regions of the sky are least or most suscepti¬ 
ble to such biases. Only after gridding many satellite 
passes together does a fuller picture emerge. Note, 
though, that given the brightness of the ORBCOMM 
satellites and the averaging discussed in 3.2, we in¬ 
terpret the statistical fluctuations as a combination 
of multipath reflections (different at the two antenna 
locations) and polarization mismatch, not as receiver 
noise. 

To characterize the behavior of these fluctuations 
as they manifest in power patterns, we combine 296 
satellite passes in the null experiment configuration 
recorded over 32 hours (spread over 4 days) into a 
measured beampattern of the reference antenna (Fig¬ 
ure 4). The measured reference antenna beampat¬ 
tern is consistent with our numerical model within 
few percent statistical scatter within 20 degrees of 
zenith, and shows modest systematic trends at the 
10% level farther out suggestive of a few degree rota¬ 
tional misalignment. This level of agreement sets an 
upper limit on beam measurement systematics due to 
environmental effects and instrumental non-idealities 
as discussed above. We thus interpret these results 
as measurements of the accuracy and precision of our 
beam measurement system in its current configura¬ 
tion. In Sec. 5, we discuss approaches to mitigating 
these systematics further. 

As a check on the reliability of our background 
estimation (used only to identify and avoid times 
of significant background power relative to satel¬ 
lite power), we plot in Figure 5 the observed back¬ 
ground level as a function of time versus that pre¬ 
dicted from the Global Sky Model (de Oliveira Costa 
2008) and our model reference antenna beampat¬ 
tern. Even neglecting the sun, the observed back¬ 
ground estimates agree qualitatively with the GSM 
at the ±0.5 dB level, close to its stated accuracy of 
±10%. Note that slight phase and amplitude dis¬ 
agreements at this level are expected as the GSM is 
generated through interpolation between sky maps 
at other frequencies, and errors are thus correlated 
over large scales. That the observed background lev¬ 
els are roughly consistent with the predicted galaxy 
power demonstrates that the ORBCOMM satellites 
are spaced sufficiently far apart in their orbits and 
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frequency bands so as to not overlap in time, which 
is crucial for our experiment. 

4. MWA Antenna Tile 


A22{0, </) + 7r/2) (see Appendix A), 


|Aii(0, (j>)\'^ (X sin^ (khcosO) 


16 

^i{k-Xi+rii 
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Having quantified the accuracy and precision of 
our setup, we proceed to a study of the beampattern 
of an MWA antenna tile (hereafter MWA tile). The 
MWA tile consists of a 4x4 grid of dual-polarization 
bowtie dipoles whose signals are combined in a delay¬ 
line beamformer. The dipoles are vertical bowties 
with dimensions 74 cm across and 84 cm on the 
diagonal, optimized to have a broad frequency re¬ 
sponse in the 100-200 MHz band. They are mounted 
on a 5 m X 5 m ground screen attached to a lev¬ 
eled wooden frame approximately 20 cm above the 
ground. The center-to-center dipole spacing is 1.1m 
and the center-to-ground-screen distance is 30 cm. 
For this experiment, we construct the ground screen 
out of 5 pieces of wire mesh (19 gauge, 0.5” spacing) 
which are crimped together every 5 inches to form a 
constant potential surface connected to earth ground. 
Each dipole has a 20 dB LNA with integrated balun. 
The beamformer sums 16 dipole signals for each of 
the two polarizations (NS and EW), producing a 
beam with Full-Width-at-Half-Maximum (FWHM) 
~ 23° at 137 MHz, with sidelobes at the —20 dB 
level. By digitally engaging delay lines in increments 
of ~450 ps to each signal pathway, the beam may 
be steered far from zenith. The delays are engaged 
through simple digital control with no amplitude or 
phase calibration needed. Slight deviation from per¬ 
fect gain and delay matching (see Sec. 4.3) across 
the beamformer pathways is one of the mismodeling 
effects this work will probe. Lastly, we note that 
beamformer adds 30 dB of gain to the summed sig¬ 
nal, to which we add 36 dB of attenuation to avoid 
saturating the ADC in our receiver. 

4.1. Model Beampatterns 

As 137 MHz is well below the half wavelength fre¬ 
quencies of the characteristic lengths of the MWA 
Dipole (202 MHz and 178 MHz), the Hertzian dipole 
model is expected to be valid, though deviations 
near the sidelobes would not be unexpected. The 
phased array and ground screen factors of the MWA 
tile are encapsulated in the array factor Aii(0, (j)) — 


Here h = 0.3 m is the height of the dipole mid¬ 
points off the ground screen, k is the direction of 
the satellite, r]i is the phase delay applied by the 
beamformer, and T) is the position of antenna i on 
the grid in the xy plane. The phase delay can be 
expressed as rn = 27r/(iix(435 ps) where di is an 
integer between 0 and 31 which is specified when 
controlling the beamformer. Lastly, we use coordi¬ 
nates where x points towards the East, y points to¬ 
wards the North, and z points towards Zenith, and 
k = ^ (sin 9 sin (j)x + sin 6 cos (j)y + cos 9z). Combin¬ 
ing Equations 2, A3, and A8 gives our analytic MWA 
tile model. 

We also compare our measurements to the more 
advanced beam model presented by Sutinjo et al. 
[2015] which includes dipole-dipole coupling effects 
and a numerical dipole model on realistically mod¬ 
eled soil. We use their average embedded element 
model, as work on a full electromagnetic coupling 
model is ongoing. Below ~ 180 MHz, the correc¬ 
tions due to dipole coupling effects are small within 
the main lobe, but potentially observable in beam 
measurement extending into the sidelobes, like ours. 

4.2. Beampattern Measurements 

Eigure 6 shows our measured MWA tile beampat¬ 
tern (top panel) when pointed towards zenith, con¬ 
structed from ~400 satellite passes recorded over ~4 
days. The beam is plotted on the HEALPIX grid dis¬ 
cussed in Sec. 3.2. We also plot the measured beam- 
pattern and our analytic model on slices through 
these polar plots on the E and H (red and blue, re¬ 
spectively) antenna planes (middle panel), as well 
as the ratio of measured over analytic model beams 
(bottom panel). The numerical model of Sutinjo 
et al. [2015] is plotted relative to the analytic model 
in the bottom panel for comparison (dashed lines), 
and should align with the data points if it explains 
the observed deviations. While measured beampat¬ 
terns are often compared with models in simple beam 
sensitivity difference plots, we view ratio plots (i.e., 
differences of dB quantities) as more relevant given 
that primary beam sky weighting during both cali¬ 
bration and primary beam correction are multiplica¬ 
tive operations. Ratio plots also highlight the off- 
zenith regions where beams are typically most poorly 
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modeled, the regions where foregrounds are most at 
risk of affecting EOR science, as discussed in Sec. 1. 

Within the half power point (~12 degrees away 
from zenith) we observe agreement with our analytic 
model beam pattern within ~5% statistical fluctua¬ 
tions. Beyond that zenith angle, towards the edges 
of the main lobe and in the few percent sidelobes, we 
observe systematic deviations away from the model 
at the dB level. We discuss these deviations and 
their patterns further in Sec. 4.3. Both positive and 
negative trends are observed. 

We have also measured the MWA tile power pat¬ 
tern at several off-zenith pointings, two of which are 
shown in Figures 7 (20 degrees East) and 8 (20 de¬ 
grees West). For these pointings, the direction of 
boresight is not in the E and H antenna planes, so 
we instead rotate the antenna E and H planes to the 
boresight direction, yielding two orthogonal planes 
crossing the off-zenith main lobe. The level of de¬ 
viation away from the model power patterns here is 
comparable with the Zenith pointing. 

4.3. Error Analysis 

A detailed analysis of the causes and effects of 
beamforming errors in MWA tiles will be presented 
in a separate paper (Neben et al., in prep). In partic¬ 
ular, that work is concerned with the finite precision 
of complex gain matching across the 16 beamformer 
signal paths (for each instrumental polarization) as 
well as tile rotation/tilt errors and dipole position 
errors. A budget of relevant systematics is estab¬ 
lished through laboratory measurements and com¬ 
pilation of manufacturer specifications and Monte 
Carlo simulations are run to propagate component 
uncertainties into direction-dependent beam power 
pattern uncertainties. Beam deviations at the level 
of 10-20% are predicted in the sidelobes and near the 
edge of the main lobe, and significantly larger in the 
nulls. 

There are, however, several sources of error in the 
beam measurements presented in this work which 
are peculiar to our MWA tile at Green Bank. Our 
ground screen is elevated 20 cm off the ground on 
a wooden frame whereas the analytic and the cross¬ 
coupling models assume it is set on the ground. This 
is expected to affect the beam pattern at low ele¬ 
vations in particular, though in a symmetric man¬ 
ner. Non-coplanarity due to ground screen sag be¬ 
tween its frame supports also complicates the ground 
screen term in Eqn. 2, and may also introduce rela¬ 


tive dipole tilts. Our ground screen is formed out of 
five rectangles of wire mesh, which are crimped to¬ 
gether every 5”. In contrast, the ground screen used 
in the deployed MWA tiles at the Murchison Radio 
Observatory in Western Australia are overlapped and 
welded, fixing a constant potential surface. 

Lastly slight tilts and rotations of our MWA tile 
due to our ±1.5° alignment precision are larger than 
those affecting deployed MWA tiles. Such alignment 
errors are most prominent where the beam changes 
rapidly with angle as it does near the edges of the 
main lobe and in the sidelobes, and are difficult 
to correct for in our beam mapping experiment as 
they upset the polarization matching with the refer¬ 
ence antenna. Numerical experiments suggest such 
alignment errors contribute systematics at the ±20% 
level. 

Many of these errors will break the ideal symme¬ 
tries of the tile beampatterns by introducing distor¬ 
tions and tilts of the main lobe and sidelobes, not 
unlike the patterns observed in the measured vs. 
model beam plots in Figures 6, 7, 8. In particu¬ 
lar, a slight main lobe widening is observed in the 
EW beams and a slight ~ 0.5° tilt is observed in the 
NS beams. These discrepancies are observed across 
all three pointings, suggesting they are due to some 
combination of the tile non-idealities discussed above 
as opposed to per-pointing gain and delay errors in 
the beamformer. 

5. Discussion 

We have used the ORBCOMM satellite constel¬ 
lation to test of the feasibility of a sky transmitter- 
based beam measurement system for low frequency 
radio interferometers. Our system compares the 
power received by an AUT to that from a well- 
modeled reference dipole, whose ratio gives the rela¬ 
tive beampattern of the AUT. The ~ 30 ORBCOMM 
satellites provide sky coverage over two thirds of the 
visible sky from our Green Bank site in a single day, 
in part due to their low earth orbits and quick orbital 
precession. Their bright signals probe deep into an¬ 
tenna sidelobes, yielding measurements over 30 dB 
of MWA tile dynamic range, even after rejecting all 
data within 20 dB of the galactic noise background. 
This is an order of magnitude improvement in beam 
measurement depth over recent source-based beam 
measurements [Colegate et al., 2014]. 

We find through our null experiment that we are 
limited by 5% statistical scatter within 20° of zenith. 



NEBEN ET AL.: MEASURING MWA BEAMPATTERNS WITH ORB COMM SATELLITES 


9 


and 10% systematics farther out. We hope to defini¬ 
tively identify these fluctuations as multipath scat¬ 
tering in the future using multi-frequency probe sig¬ 
nals to investigate their frequency dependence. The 
time scale of multipath fluctuations is set by the 
satellite’s motion through the frequency-dependent 
interference pattern set up on the ground by nearby 
reflecting structures. Finer antenna alignment at the 
sub-degree level will also help mitigate systematics 
near the edges of the main lobe. 

We have used this prototype system to conduct 
the first measurements of an MWA tile beampattern 
over a large field of view including the main lobe and 
primary sidelobes. We find good agreement with a 
simple array factor-based short dipole model within 
the tile FWHM (23 degrees across at 137 MHz), but 
observe ~ dB level systematic deviations towards the 
edge of the main lobe {9 ~ 20°) and in the sidelobes. 
These deviations are larger than the 10% systematics 
observed in our null experiment, and in principle rep¬ 
resent the beam modeling errors we originally sought 
to measure. 

However, several considerations prevent our in¬ 
terpretation of these deviations as inherent in the 
MWA tile design, and thus representative of the 128 
deployed antennas at the Murchison Radio Obser¬ 
vatory in Western Australia. Sutinjo et al. [2015] 
show MWA imaging results which are consistent with 
an advanced MWA tile model including mutual cou¬ 
pling, however our measured beampatterns appear 
no more consistent with this model than with our 
simple analytic one. Indeed the deviations we ob¬ 
serve lack the symmetry expected from a beam mod¬ 
eling error due only to insufficiently precise modeling 
of the sort the cross coupling model is designed to 
correct. In that case, the zenith pointing EW and 
NS beampatterns should coincide after a 90° rota¬ 
tion, and each should exhibit 180° rotational sym¬ 
metry. Additionally, as discussed in Sec. 4.3, known 
MWA beamforming errors are predicted to be at 
most ~ 20% towards the edge of the main lobe and in 
sidelobes. However, there exist additional errors pe¬ 
culiar to our MWA tile set up in Green Bank which 
are more difficult to quantify (e.g., tile and dipole 
tilts, imperfect ground, ground screen sag, electric 
potential non-constancy). We thus view this work 
less as a measurement of the MWA tile beampattern 
and more as a demonstration of the power of the 
ORBCOMM technique for identifying beam model¬ 
ing errors (i.e., deviations of an as-built antenna from 


from its ideal model) down to the -30 dB level of 
the beam. We expect that further tests using probe 
signals from drones or multi-frequency satellites will 
both hone our understanding of the technique and 
set tighter constraints on beam models of the MWA 
tile. 

Still unknown, though, are the effects of these er¬ 
rors on ongoing MWA Epoch of Reionization power 
spectrum measurements, which will depend both on 
their frequency dependence and the degree to which 
they limit calibration fidelity. This work has begun 
to probe these effects in a way that imaging can¬ 
not, as beamforming errors tend to average out when 
forming an image with many MWA tiles. However, 
such averaging of beamforming errors will be less per¬ 
fect when forming sky power spectra because differ¬ 
ent antennas probe different regions of the uv plane, 
depending on which baselines they are part of. In¬ 
deed as proposed by Morales and Matejek [2009], in¬ 
terferometric imaging algorithms taking per-antenna 
primary beams into account [Sullivan et al., 2012; 
Dillon et al., 2014b] may be critical in order to ac¬ 
cess the Epoch of Reionization. 

Appendix A: Measurement of the 
unpolarized beampattern 

In general the response of an antenna to unpo¬ 
larized radiation (e.g., thermal emission) is different 
from its response to polarized radiation (e.g., satellite 
signals). We show how the unpolarized beampattern 
(defined below) of the AUT may be measured despite 
any polarization of the satellite probe signal, assum¬ 
ing that both the AUT and the reference dipole have 
the same polarization response. 

In general, the voltage responses of these antennas 
to radiation from (9, (p) relate to the two incident sky 
polarizations as 

Uaut = are (Al) 

Uref = RE. (A2) 

where V = (y^) and E = (eI)- We use coordinates 
where x points to the East, and y points to the North. 
Spherical unit vectors 9 and (j) point in the directions 
of increasing 9 and </>, respectively. These unit vec¬ 
tors are thus functions of those angles, though here 
we consider radiation incident only from the single di¬ 
rection {9, (f). Further, R is a matrix which converts 
from sky polarization to instrument polarization and 
drops the z component (it is a subset of a rotation 
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matrix), 

cos0sine/) coscj) \ , 

cos 9 cos (j) — sin </> / ^ 

Assuming both antennas have the same polariza¬ 
tion response, in the sense that the x (y) oriented 
dipoles respond only to x (y) polarized radiation, A 
is a diagonal matrix which is a function only of angle 
on the sky and does not mix instrument polariza¬ 
tions. The physical origin of A is the array factor 
of the MWA tile, as well as effects of ground screen, 
antenna geometry, and any dipole cross-coupling. In 
the ideal case, Aii{9, (j)) = A22(9, (j) + 7r/2), though 
our measurement does not assume this. 

Consider the x oriented antennas as an example. 
The received powers in response to polarized radia¬ 
tion are 

Paut.x = |Ann[A’?i(|A,|2) + 2i?ni?i2W^^0*) 

+Rl2m?)] (A4) 

Pref.x = Rl^{\E0\^) +2R^iRi2^{EeE^*) 

+Rl2{\E^?)- (A5) 

However, if the incident radiation is unpolarized 
with intensity I, then Eg and E^f, are uncorrelated 
and equal in power and may be added together 
as powers to give the total received power (/ = 
2{\Eg\'^) = 2(|i?0p)). The received powers are then 


Paut,x = |Ann(i?2^+i?2^)//2 (A6) 

P.ei,. = {R\^ + Rl2)Il2, (A7) 

and so the unpolarized x beampatterns are given by 
Saut.x = |Ann(i?2^+ii-2^)/2 (A8) 

Bref.x = (i??! + i??2)/2. (A9) 


To measure Haut.x; we thus need only An, which 
can be computed from the ratio of Paut,x and Pref.x, 
as R is known (see Eqns. A4 and A5). This is exactly 
the procedure outlined in Sec. 3.1. 

Acknowledgments. This work was supported by 
NSF grant AST-0821321, the Marble Astrophysics Fund, 
and the MIT School of Science. We thank Pat Klima, 
Bang Dinh Nhan, and the staff at the National Radio 
Astronomy Observatory-Green Bank for assistance in 
setting up and debugging our experiment, and Haox- 
uan Zheng, Aaron Ewall-Wice, Josh Dillon, Lu Feng, 
and Daniel Jacobs for helpful discussions. We thank 
Nithyanandan Thyagarajan, Josh Dillon, Adrian Sutinjo, 
Tim Colegate, and the anonymous referees for very help¬ 
ful comments on our manuscript. 

This scientific work makes use of the Murchison Radio- 
astronomy Observatory, operated by CSIRO. We ac¬ 
knowledge the Wajarri Yamatji people as the traditional 


owners of the Observatory site. Support for the MWA 
comes from the U.S. National Science Foundation (grants 
AST-0457585, PHY-0835713, CAREER-0847753, and 
AST-0908884), the Australian Research Council (LIEF 
grants LE0775621 and LE0882938), the U.S. Air Force 
Office of Scientific Research (grant FA9550-0510247), and 
the Centre for All-sky Astrophysics (an Australian Re¬ 
search Council Centre of Excellence funded by grant 
CE110001020). Support is also provided by the Smithso¬ 
nian Astrophysical Observatory, the MIT School of Sci¬ 
ence, the Raman Research Institute, the Australian Na¬ 
tional University, and the Victoria University of Welling¬ 
ton (via grant MED-E1799 from the New Zealand Min¬ 
istry of Economic Development and an IBM Shared Uni¬ 
versity Research Grant). The Australian Federal gov¬ 
ernment provides additional support via the Common¬ 
wealth Scientific and Industrial Research Organisation 
(CSIRO), National Collaborative Research Infrastructure 
Strategy, Education Investment Fund, and the Australia 
India Strategic Research Fund, and Astronomy Australia 
Limited, under contract to Curtin University. We ac¬ 
knowledge the iVEC Petabyte Data Store, the Initiative 
in Innovative Computing and the CUDA Center for Ex¬ 
cellence sponsored by NVIDIA at Harvard University, 
and the International Centre for Radio Astronomy Re¬ 
search (ICRAR), a Joint Venture of Curtin University 
and The University of Western Australia, funded by the 
Western Australian State government. Data on which fig¬ 
ures and tables herein are based may be obtained by con¬ 
tacting the corresponding author Abraham Neben (abra- 
hamn@mit.edu). 

Notes 

1. http://www.qsl.net/kd2bd/predict.html 

2. http://www.celestrak.com/NORAD/elements/orbcomm.txt 

3. https://www.cst.com/Products/CSTMWS 


References 

Asad, K. M. B., L. V. E. Koopmans, V. Jelic, V. N. 
Pandey, A. Ghosh, F. B. Abdalla, G. Bernardi, M. A. 
Brentjens, A. G. de Bruyn, S. Bus, B. Ciardi, E. Chap¬ 
man, S. Daiboo, E. R. Fernandez, G. Harker, 1. T. 
Iliev, H. Jensen, O. Martinez-Rubi, G. Mellema, 
M. Mevius, A. R. Offringa, A. H. Patil, J. Schaye, 

R. M. Thomas, S. van der Tol, H. K. Vedantham, 

S. Yatawatta, and S. Zaroubi (2015), Polarization leak¬ 
age in Epoch of Reionization windows: 1. LOFAR ob¬ 
servations of the 3C196 field, ArXiv e-prints. 

Beardsley, A. P., B. J. Hazelton, M. F. Morales, W. Ar¬ 
cus, D. Barnes, G. Bernardi, J. D. Bowman, F. H. 
Briggs, J. D. Bunton, R. J. Cappallo, B. E. Corey, 

A. Deshpande, L. deSouza, D. Emrich, B. M. Gaensler, 
R. Goeke, L. J. Greenhill, D. Herne, J. N. Hewitt, 
M. Johnston-Hollitt, D. L. Kaplan, J. C. Kasper, 

B. B. Kincaid, R. Koenig, E. Kratzenberg, C. J. Lons¬ 
dale, M. J. Lynch, S. R. MeWhirter, D. A. Mitchell, 
E. Morgan, D. Oberoi, S. M. Ord, J. Pathikulangara, 

T. Prabu, R. A. Remillard, A. E. E. Rogers, A. Roshi, 



NEBEN ET AL.: MEASURING MWA BEAMPATTERNS WITH ORB COMM SATELLITES 


11 


J. E. Salah, R. J. Sault, S. N. Udaya, K. S. Srivani, 
J. Stevens, R. Subrahmanyan, S. J. Tingay, R. B. 
Wayth, M. Waterson, R. L. Webster, A. R. Whit¬ 
ney, A. Williams, C. L. Williams, and J. S. B. Wyithe 
(2013), The EoR sensitivity of the Murchison Wide- 
held Array, Monthly Notices of The Royal Astronomi¬ 
cal Society, 429, L5-L9, doi:10.1093/mnrasl/sls013. 
Bernard!, G., L. J. Greenhill, D. A. Mitchell, S. M. Ord, 

B. J. Hazelton, B. M. Gaensler, A. de Oliveira-Costa, 
M. E. Morales, N. Udaya Shankar, R. Subrahmanyan, 
R. B. Wayth, E. Lenc, C. L. Williams, W. Arcus, B. S. 
Arora, D. G. Barnes, J. D. Bowman, F. H. Briggs, 
J. D. Bunton, R. J. Cappallo, B. E. Corey, A. Desh- 
pande, L. deSouza, D. Emrich, R. Goeke, D. Herne, 
J. N. Hewitt, M. Johnston-Hollitt, D. Kaplan, J. C. 
Kasper, B. B. Kincaid, R. Koenig, E. Kratzenberg, 

C. J. Lonsdale, M. J. Lynch, S. R. McWhirter, E. Mor¬ 
gan, D. Oberoi, J. Pathikulangara, T. Prabu, R. A. 
Remillard, A. E. E. Rogers, A. Roshi, J. E. Salah, R. J. 
Sault, K. S. Srivani, J. Stevens, S. J. Tingay, M. Wa¬ 
terson, R. L. Webster, A. R. Whitney, A. Williams, 
and J. S. B. Wyithe (2013), A 189 MHz, 2400 deg^ 
Polarization Survey with the Murchison Wideheld Ar¬ 
ray 32-element Prototype, Astrophys. J., , 111, 105, 
doi:10.1088/0004-637X/771/2/105. 

Bowman, J. D., D. G. Barnes, E. H. Briggs, B. E. Corey, 
M. J. Lynch, N. D. R. Bhat, R. J. Cappallo, S. S. 
Doeleman, B. J. Fanous, D. Herne, J. N. Hewitt, 
C. Johnston, J. C. Kasper, J. Kocz, E. Kratzenberg, 

C. J. Lonsdale, M. F. Morales, D. Oberoi, J. E. Salah, 

B. Stansby, J. Stevens, G. Torr, R. Wayth, R. L. Web¬ 
ster, and J. S. B. Wyithe (2007), Field deployment of 
prototype antenna tiles for the mileura wideheld array 
low frequency demonstrator. The Astronomical Jour¬ 
nal, 133(4:), 1505. 

Bowman, J. D., 1. Cairns, D. L. Kaplan, T. Murphy, 

D. Oberoi, L. Staveley-Smith, W. Arcus, D. G. Barnes, 
G. Bernard!, F. H. Briggs, S. Brown, J. D. Bunton, 
A. J. Burgasser, R. J. Cappallo, S. Chatterjee, B. E. 
Corey, A. Coster, A. Deshpande, L. deSouza, D. Em- 
rich, P. Erickson, R. F. Goeke, B. M. Gaensler, L. J. 
Greenhill, L. Harvey-Smith, B. J. Hazelton, D. Herne, 
J. N. Hewitt, M. Johnston-Hollitt, J. C. Kasper, B. B. 
Kincaid, R. Koenig, E. Kratzenberg, C. J. Lonsdale, 
M. J. Lynch, L. D. Matthews, S. R. McWhirter, D. A. 
Mitchell, M. E. Morales, E. H. Morgan, S. M. Ord, 
J. Pathikulangara, T. Prabu, R. A. Remillard, T. Ro- 
bishaw, A. E. E. Rogers, A. A. Roshi, J. E. Salah, 
R. J. Sault, N. U. Shankar, K. S. Srivani, J. B. Stevens, 
R. Subrahmanyan, S. J. Tingay, R. B. Wayth, M. Wa¬ 
terson, R. L. Webster, A. R. Whitney, A. J. Williams, 

C. L. Williams, and J. S. B. Wyithe (2013), Science 
with the murchison wideheld array, TrottObserving- 
Modes, 30, doi:10.1017/pas.2013.009. 

Bradley, R., and P. Ries (2008), Mapping the beam of 
a low frequency antenna using satellite signals (con¬ 
ference presentation, invited), national Radio Science 
Meeting, Boulder, CO. 

Brueckmann, H. (1963), Antenna pattern measurement 
by satellite, IEEE Transactions on Antennas and Pro¬ 
pagation, p. 143. 


Butler, B. (1998), Options for VLBA antennan surface 
measurement, Tech, rep.. National Radio Astronomy 
Observatory, VLBA Test Memo. No. 57. 

Colegate, T., A. Sutinjo, P. Hall, S. Padhi, R. Wayth, 
J. de Vaate, B. Crosse, D. Emrich, A. Eaulkner, 
N. Hurley-Walker, E. de Lera Acedo, B. Juswardy, 
N. Razavi-Ghods, S. Tingay, and A. Williams 
(2014), Antenna array characterization via ra¬ 
dio interferometry observation of astronomical 
sources, in Antenna Measurements Applications 
(GAMA), 2014 IEEE Conference on, pp. 1-4, 
doi:10.1109/CAMA.2014.7003436. 

Czekala, L, and R. Bradley (2010a), Calibrating astro¬ 
nomical antenna arrays with man-made satellites (con¬ 
ference presentation), american Institute of Aeronau¬ 
tics and Astronautics (AIAA), Region I-MA Student 
Conference, Blacksburg, VA. 

Czekala, L, and R. Bradley (2010b), Calibrating astro¬ 
nomical antenna arrays with orbiting data satellites 
(conference presentation), atlantic Coast Conference 
(ACC) ‘Meeting of the Minds’ Conference, Atlanta, 
GA. 

Datta, A., J. D. Bowman, and C. L. Carilli (2010), Bright 
Source Subtraction Requirements for Redshifted 21 cm 
Measurements, ApJ, 724, 526-538, doi:10.1088/0004- 
637X/724/1/526. 

De Oliveira-Costa, A., M. Tegmark, B. M. Gaensler, 
J. Jonas, T. L. Landecker, and P. Reich (2008), 
A model of dihuse galactic radio emission from 10 
mhz to 100 ghz. Monthly Notices of the Royal Astro¬ 
nomical Society, 388(1), 247-260, doi:10.1111/j.l365- 
2966.2008.13376.x. 

Deguchi, H., and et al. (1993), Radio holographic metrol¬ 
ogy with best-ht panel model of the nobeyama 45-m 
telescope, lEICE Trans. Comm., 37, 1492-1499. 

Dillon, J. S., A. Liu, C. L. Williams, J. N. Hewitt, 
M. Tegmark, E. H. Morgan, A. M. Levine, M. F. 
Morales, S. J. Tingay, G. Bernard!, J. D. Bow¬ 
man, F. H. Briggs, R. C. Cappallo, D. Emrich, 
D. A. Mitchell, D. Oberoi, T. Prabu, R. Wayth, and 
R. L. Webster (2014a), Overcoming real-world obsta¬ 
cles in 21 cm power spectrum estimation: A method 
demonstration and results from early Murchison Wide- 
held Array data, Rhys. Rev. D, 89(2), 023002, doi: 
10.1103/PhysRevD.89.023002. 

Dillon, J. S., M. Tegmark, A. Liu, A. Ewall-Wice, J. N. 
Hewitt, M. F. Morales, A. R. Neben, A. R. Parsons, 
and H. Zheng (2014b), Mapmaking for Precision 21 
cm Cosmology, ArXiv e-prints. 

Fukao, S., T. Sato, and S. Kato (1985), Monitoring of 
the mu radar antenna pattern by satellite OHZORA 
(OXOS-C), J. Geomag. Geoelectr., 37, 431. 

Furlanetto, S., S. Oh, and F. Briggs (2006), Cosmology 
at low frequencies: The 21 cm transition and the high- 
redshift Universe, Physics Reports, 433, 181-301, doi: 
10.1016/j.physrep.2006.08.002. 

Godwin, and et al. (1986), Improvement of the ehelsberg 
100 meter telescope based on holographic rehector sur¬ 
face measurements, Astron. Astro., 167, 390-394. 

Gorski, K. M., E. Hivon, A. J. Banday, B. D. Wan- 
delt, F. K. Hansen, M. Reinecke, and M. Bartelmann 
(2005), Healpix: A framework for high-resolution dis¬ 
cretization and fast analysis of data distributed on the 
sphere. The Astrophysical Journal, 622(2), 759. 



12 


NEBEN ET AL.: MEASURING MWA BEAMPATTERNS WITH ORBCOMM SATELLITES 


Hall, P. J. (2005), The Square Kilometre Array: An En¬ 
gineering Perspective. 

Harp, G., R. Ackermann, Z. Nadler, S. Blair, M. Davis, 
M. Wright, J. Eorster, D. DeBoer, W. Welch, S. Atkin¬ 
son, D. Backer, P. Backus, W. Barott, A. Bauer- 
meister, L. Blitz, D.-J. Bock, G. Bower, T. Bradford, 
C. Cheng, S. Croft, M. Dexter, J. Dreher, G. Engargi- 
ola, E. Fields, C. Heiles, T. Heifer, J. Jordan, S. Jor¬ 
gensen, T. Kilsdonk, C. Gutierrez-Kraybill, G. Keat¬ 
ing, C. Law, J. Lugten, D. MacMahon, P. McMahon, 
O. Milgrome, A. Siemion, K. Smolek, D. Thornton, 
T. Pierson, K. Randall, J. Ross, S. Shostak, J. Tarter, 
L. Urry, D. Werthimer, P. Williams, and D. Whysong 
(2011), Primary beam and dish surface characteriza¬ 
tion at the alien telescope array by radio hologra¬ 
phy, Antennas and Propagation, IEEE Transactions 
on, 59(6), 2004-2021, doi:10.1109/TAP.2011.2122214. 

Hurtado, M., H. Lorente, and C. Muravchik (2001), An 
in-situ method to measure the radiation pattern of a 
gps receiving antenna, IEEE Transactions on Instru¬ 
mentation and Measurement, 50, 846. 

Jacobs, D., J. Bowman, and J. Aguirre (2013), The pre¬ 
cision and accuracy of early epoch of reionization fore¬ 
ground models: comparing MWA and PAPER 32- 
antenna source catalogs. The Astrophysical Journal, 
769(1), 5. 

Jelic, V., S. Zaroubi, P. Labropoulos, G. Bernard!, 

A. G. de Bruyn, and L. V. E. Koopmans (2010), 
Realistic simulations of the galactic polarized fore¬ 
ground: consequences for 21-cm reionization detection 
experiments. Monthly Notices of the Royal Astronom¬ 
ical Society, 409(4), 1647-1659, doi:10.1111/j.l365- 
2966.2010.17407.x. 

Lasenby, A. (1985), Surface measurements of the NRAO 
12m antenna, in Proc. ESO-IRAM-Onsala Workshop 
on ’(Sub)Milimeter Astronomy’, edited by P. Shaver 
and K. Kjar, pp. 77-83. 

Law, D. C., J. R. Khorrami, W. B. Sessions, and M. K. 
Shanahan (1997), Radiation patterns of a large UHE 
phased-array antenna: A comparison of measurements 
using satellite repeaters and patterns derived from 
measurements of antenna current distributions, IEEE 
Transactions on Antennas and Propagation Magazine, 
39, 88. 

Liu, A., A. R. Parsons, and C. M. Trott (2014a), 
Epoch of reionization window. 1. Mathematical 
formalism, Phys. Rev. D, 90(2), 023018, doi: 

10.1103/PhysRevD.90.023018. 

Liu, A., A. R. Parsons, and C. M. Trott (2014b), Epoch 
of reionization window. H. Statistical methods for fore¬ 
ground wedge reduction, Phys. Rev. D, 90(2), 023019, 
doi:10.1103/PhysRevD.90.023019. 

Loeb, A., and S. R. Furlanetto (2013), The First Galaxies 
In The Universe, Princeton University Press, Prince¬ 
ton, NJ. 

Lonsdale, C., R. Cappallo, M. Morales, F. Briggs, 
L. Benkevitch, J. Bowman, J. Bunton, S. Burns, 

B. Gorey, L. deSouza, S. Doeleman, M. Derome, 
A. Deshpande, M. Gopala, L. Greenhill, D. Herne, 
J. Hewitt, P. Kamini, J. Kasper, B. Kincaid, 
J. Kocz, E. Kowald, E. Kratzenberg, D. Kumar, 


M. Lynch, S. Madhavi, M. Matejek, D. Mitchell, 
E. Morgan, D. Oberoi, S. Ord, J. Pathikulangara, 
T. Prabu, A. Rogers, A. Roshi, J. Salah, R. Sault, 

N. Shankar, K. Srivani, J. Stevens, S. Tingay, 
A. Vaccarella, M. Waterson, R. Wayth, R. Webster, 
A. Whitney, A. Williams, and G. Williams (2009), 
The murchison widefield array: Design overview. 
Proceedings of the IEEE, 97(8), 1497-1506, doi: 
10.1109/JPROG.2009.2017564. 

Moore, D., J. E. Aguirre, A. Parsons, Z. Ali, R. Bradley, 
G. Carilli, D. DeBoer, M. Dexter, N. Gugliucci, D. Ja¬ 
cobs, P. Klima, A. Liu, D. MacMahon, J. Manley, 
J. Pober, 1. Stefan, and W. Walbrugh (2015), New 
Limits on Polarized Power Spectra at 126 and 164 
MHz: Relevance to Epoch of Reionization Measure¬ 
ments, ArXiv e-prints. 

Moore, D. F., J. E. Aguirre, A. R. Parsons, D. C. Ja¬ 
cobs, and J. C. Pober (2013), The effects of polar¬ 
ized foregrounds on 21cm epoch of reionization power 
spectrum measurements, The Astrophysical Journal, 
769(2), 154. 

Morales, M. E., and M. Matejek (2009), Software holog¬ 
raphy: interferometric data analysis for the chal¬ 
lenges of next generation observatories, Mon. Not. R. 
Astron. Soc., , fOO, 1814-1820, doi:10.1111/j.l365- 
2966.2009.15537.x. 

Morales, M. E., and J. S. B. Wyithe (2010), Reioniza¬ 
tion and Cosmology with 21-cm Fluctuations, Annual 
Reviews of Astronomy and Astrophysics, 48, 127-171, 
doi: 10.1146/annurev-ast ro-081309-130936. 

Morales, M. F., B. Hazelton, 1. Sullivan, and A. Beard¬ 
sley (2012), Four Fundamental Foreground Power 
Spectrum Shapes for 21 cm Cosmology Observa¬ 
tions, Astrophys. J., , 752, 137, doi:10.1088/0004- 
637X/752/2/137. 

Neben, A., R. Bradley, and J. Hewitt (2014), Measur¬ 
ing antenna tile beams with high dynamic range for 
the murchison widefield array (MWA) using 137 mhz 
ORBCOMM satellites, american Astronomical Society 
Meeting, Boston, MA. 

Parsons, A., J. Pober, M. McQuinn, D. Jacobs, 
and J. Aguirre (2012), A Sensitivity and Array- 
configuration Study for Measuring the Power Spec¬ 
trum of 21 cm Emission from Reionization, Astrophys. 
J., , 753, 81, doi:10.1088/0004-637X/753/l/81. 

Parsons, A. R., D. C. Backer, G. S. Foster, M. C. H. 
Wright, R. F. Bradley, N. E. Gugliucci, C. R. 
Parashare, E. E. Benoit, J. E. Aguirre, D. C. Jacobs, 

C. L. Carilli, D. Herne, M. J. Lynch, J. R. Manley, and 

D. J. Werthimer (2010), The precision array for prob¬ 
ing the epoch of re-ionization: Eight station results. 
The Astronomical Journal, 139(4), 1468. 

Parsons, A. R., A. Liu, J. E. Aguirre, Z. S. Ali, R. F. 
Bradley, et al. (2014), New Limits on 21cm EoR From 
PAPER-32 Consistent with an X-Ray Heated IGM 
at z=7.7, Astrophys.J., 788, 106, doi:10.1088/0004- 
637X/788/2/106. 

Pober, J., A. Parsons, D. Jacobs, J. Aguirra, R. Bradley, 
C. Carilli, N. Gugliucci, D. Moore, and C. Parashare 
(2012), A technique for primary beam calibration of 
drift-scanning, wide-held antenna elements, The As¬ 
tronomical Journal, 143(2), 53. 



NEBEN ET AL.: MEASURING MWA BEAMPATTERNS WITH ORB COMM SATELLITES 


13 


Pober, J. C., A. R. Parsons, J. E. Aguirre, Z. Ali, R. F. 
Bradley, C. L. Carilli, D. DeBoer, M. Dexter, N. E. 
Gugliucci, D. C. Jacobs, D. MacMahon, J. Manley, 

D. F. Moore, 1. 1. Stefan, and W. P. Walbrugh (2013), 
Opening the 21cni EoR Window: Measurements of 
Foreground Isolation with PAPER, ArXiv e-prints. 

Pober, J. C., A. Lin, J. S. Dillon, J. E. Agnirre, J. D. Bow¬ 
man, R. F. Bradley, C. L. Carilli, D. R. DeBoer, J. N. 
Hewitt, D. C. Jacobs, M. McQuinn, M. F. Morales, 

A. R. Parsons, M. Tegmark, and D. J. Werthimer 
(2014), What next-generation 21 cm power spectrnm 
measurements can teach ns about the epoch of reion¬ 
ization, The Astrophysical Journal, 782(2), 66. 
Pritchard, J. R., and A. Loeb (2012), 21 cm cosmology 
in the 21st century. Reports on Progress in Physics, 
75(8), 086,901, doi:10.1088/0034-4885/75/8/086901. 
Ries, P. (2007), Beam patterns of low freqnnecy radio 
antennas using LEO satellites, M.S. thesis. University 
of Virginia, Charlottesville. 

Rochblatt, D., and B. Seidel (1992), Microwave antenna 
holography, IEEE Trans. Mic. Theory Tech., MTT-40, 
1294-1300. 

Sullivan, 1. S., M. F. Morales, B. J. Hazelton, W. Ar¬ 
ens, D. Barnes, G. Bernard!, F. H. Briggs, J. D. Bow¬ 
man, J. D. Bunton, R. J. Cappallo, B. E. Corey, 

A. Deshpande, L. deSouza, D. Emrich, B. M. Gaensler, 

R. Goeke, L. J. Greenhill, D. Herne, J. N. Hewitt, 
M. Johnston-Hollitt, D. L. Kaplan, J. C. Kasper, 

B. B. Kincaid, R. Koenig, E. Kratzenberg, C. J. Lons¬ 
dale, M. J. Lynch, S. R. MeWhirter, D. A. Mitchell, 

E. Morgan, D. Oberoi, S. M. Ord, J. Pathikulangara, 
T. Prabu, R. A. Remillard, A. E. E. Rogers, A. Roshi, 
J. E. Salah, R. J. Sault, N. Udaya Shankar, K. S. 
Srivani, J. Stevens, R. Subrahmanyan, S. J. Tingay, 

R. B. Wayth, M. Waterson, R. L. Webster, A. R. 
Whitney, A. Williams, C. L. Williams, and J. S. B. 
Wyithe (2012), Fast Holographic Deconvolution: A 
New Techniqne for Precision Radio Interferometry, 
The Astrophysical Journal, 759, 17, doi:10.1088/0004- 
637X/759/1/17. 

Sutinjo, A., J. O’Sullivan, E. Lenc, R. B. Wayth, 

S. Padhi, P. Hall, and S. J. Tingay (2015), Understand¬ 
ing instrumental stokes leakage in murchison widefield 
array polarimetry. Radio Science, 55(1), 52-65, doi: 
10.1002/2014RS005517. 

Thyagarajan, N., D. J. Helfand, R. L. White, and R. H. 
Becker (2011), Variable and transient radio sources in 
the FIRST survey. The Astrophysical Journal, 7/^(1), 
49. 

Thyagarajan, N., N. U. Shankar, R. Subrahmanyan, 
W. Arens, G. Bernard!, J. D. Bowman, F. Briggs, 
J. D. Bunton, R. J. Cappallo, B. E. Corey, L. deS¬ 
ouza, D. Emrich, B. M. Gaensler, R. F. Goeke, L. J. 
Greenhill, B. J. Hazelton, D. Herne, J. N. Hewitt, 
M. Johnston-Hollitt, D. L. Kaplan, J. C. Kasper, 
B. B. Kincaid, R. Koenig, E. Kratzenberg, C. J. Lons¬ 
dale, M. J. Lynch, S. R. MeWhirter, D. A. Mitchell, 
M. F. Morales, E. H. Morgan, D. Oberoi, S. M. Ord, 
J. Pathikulangara, R. A. Remillard, A. E. E. Rogers, 
D. A. Roshi, J. E. Salah, R. J. Sault, K. S. Srivani, 
J. B. Stevens, P. Thiagaraj, S. J. Tingay, R. B. Wayth, 
M. Waterson, R. L. Webster, A. R. Whitney, A. J. 
Williams, C. L. Williams, and J. S. B. Wyithe (2013), 
A study of fundamental limitations to statistical detec¬ 
tion of redshifted h i from the epoch of reionization. 
The Astrophysical Journal, 775(1), 6. 


Thyagarajan, N., D. C. Jacobs, J. D. Bowman, N. Barry, 

A. P. Beardsley, G. Bernard!, F. Briggs, R. J. Cap¬ 
pallo, P. Carroll, B. E. Corey, A. de Oliveira-Costa, 
J. S. Dillon, D. Emrich, A. Ewall-Wice, L. Feng, 
R. Goeke, L. J. Greenhill, B. J. Hazelton, J. N. Hewitt, 
N. Hurley-Walker, M. Johnston-Hollitt, D. L. Kaplan, 
J. C. Kasper, H.-S. Kim, P. Kittiwisit, E. Kratzenberg, 
E. Lenc, J. Line, A. Loeb, C. J. Lonsdale, M. J. Lynch, 

B. McKinley, S. R. MeWhirter, D. A. Mitchell, M. F. 
Morales, E. Morgan, A. R. Neben, D. Oberoi, A. R. 
Offringa, S. M. Ord, S. Paul, B. Pindor, J. C. Pober, 

T. Prabu, P. Procopio, J. Riding, A. E. E. Rogers, 

A. Roshi, N. U. Shankar, S. K. Sethi, K. S. Srivani, 

R. Subrahmanyan, 1. S. Sullivan, M. Tegmark, S. J. 
Tingay, C. M. Trott, M. Waterson, R. B. Wayth, R. L. 
Webster, A. R. Whitney, A. Williams, C. L. Williams, 

C. Wu, and J. S. B. Wyithe (2015), Foregrounds in 
wide-held redshifted 21 cm power spectra, The Astro- 
physical Journal, 50/(1), 14. 

Tingay, S. J., R. Goeke, J. D. Bowman, D. Emrich, 

S. M. Ord, D. A. Mitchell, M. F. Morales, T. Booler, 

B. Crosse, R. B. Wayth, C. J. Lonsdale, S. Tremblay, 

D. Pallot, T. Colegate, A. Wicenec, N. Kudryavtseva, 
W. Arcus, D. Barnes, G. Bernard!, F. Briggs, S. Burns, 
J. D. Bunton, R. J. Cappallo, B. E. Corey, A. Desh¬ 
pande, L. Desouza, B. M. Gaensler, L. J. Greenhill, 
P. J. Hall, B. J. Hazelton, D. Herne, J. N. Hewitt, 
M. Johnston-Hollitt, D. L. Kaplan, J. C. Kasper, B. B. 
Kincaid, R. Koenig, E. Kratzenberg, M. J. Lynch, 
B. Mckinley, S. R. Mewhirter, E. Morgan, D. Oberoi, 
J. Pathikulangara, T. Prabu, R. A. Remillard, A. E. E. 
Rogers, A. Roshi, J. E. Salah, R. J. Sault, N. Udaya- 
Shankar, F. Schlagenhaufer, K. S. Srivani, J. Stevens, 
R. Subrahmanyan, M. Waterson, R. L. Webster, A. R. 
Whitney, A. Williams, C. L. Williams, and J. S. B. 
Wyithe (2013), The murchison wideheld array: The 
square kilometre array precursor at low radio frequen¬ 
cies, Publications of the Astronomical Society of Aus¬ 
tralia, 30, doi:10.1017/pasa.2012.007. 

Trott, C. M., R. B. Wayth, and S. J. Tingay (2012), The 
Impact of Point-source Subtraction Residuals on 21 
cm Epoch of Reionization Estimation, Astrophys. J., 
, 757, 101, doi:10.1088/0004-637X/757/l/101. 

van Haarlem, M. R, M. W. Wise, A. W. Gunst, G. Heald, 
J. P. McKean, J. W. T. Hessels, A. G. de Bruyn, 
R. Nijboer, J. Swinbank, R. Fallows, M. Brentjens, 

A. Nelles, R. Beck, H. Falcke, R. Fender, J. Horandel, 

L. V. E. Koopmans, G. Mann, G. Miley, H. Rottgering, 

B. W. Stappers, R. A. M. J. Wijers, S. Zaroubi, 

M. van den Akker, A. Alexov, J. Anderson, K. An¬ 
derson, A. van Ardenne, M. Arts, A. Asgekar, 1. M. 
Avruch, F. Batejat, L. Bahren, M. E. Bell, M. R. 
Bell, 1. van Bemmel, P. Bennema, M. J. Bentum, 
G. Bernard!, P. Best, L. Birzan, A. Bonafede, A.-J. 
Boonstra, R. Braun, J. Bregman, F. Breitling, R. H. 
van de Brink, J. Broderick, P. C. Broekema, W. N. 
Brouw, M. Briiggen, H. R. Butcher, W. van Cap- 
pellen, B. Ciardi, T. Coenen, J. Conway, A. Coolen, 



14 


NEBEN ET AL.: MEASURING MWA BEAMPATTERNS WITH ORBCOMM SATELLITES 


A. Corstanje, S. Damstra, O. Davies, A. T. Deller, R.- 
J. Dettmar, G. van Diepen, K. Dijkstra, P. Donker, 
A. Doorduin, J. Dromer, M. Drost, A. van Duin, 
J. Eisloffel, J. van Enst, C. Ferrari, W. Frieswijk, 
H. Gankema, M. A. Garrett, F. de Gasperin, M. Ger- 
bers, E. de Gens, J.-M. Griefimeier, T. Grit, P. Grup- 
pen, J. P. Hamaker, T. Hassall, M. Hoeft, H. A. 
Holties, A. Horneffer, A. van der Horst, A. van 
Houwelingen, A. Hnijgen, M. lacobelli, H. Interna, 
N. Jackson, V. Jelic, A. de Jong, E. Juette, D. Kant, 
A. Karastergiou, A. Koers, H. Kollen, V. 1. Kon¬ 
dratiev, E. Kooistra, Y. Koopman, A. Koster, M. Ku- 
niyoshi, M. Kramer, G. Kuper, P. Lambropoulos, 
C. Law, J. van Leenwen, J. Lemaitre, M. Loose, 
P. Maat, G. Macario, S. Markoff, J. Masters, R. A. Mc- 
Fadden, D. McKay-Bukowski, H. Meijering, H. Menl- 
man, M. Mevius, E. Middelberg, R. Millenaar, J. G. A. 
Miller-Jones, R. N. Mohan, J. D. Mol, J. Morawietz, 

R. Morganti, D. D. Mulcahy, E. Mulder, H. Munk, 

L. Nieuwenhuis, R. van Nienwpoort, J. E. Noordam, 

M. Norden, A. Noutsos, A. R. Offringa, H. Olof- 
sson, A. Omar, E. Orrir, R. Overeem, H. Paas, 
M. Pandey-Pommier, V. N. Pandey, R. Pizzo, A. Po- 
latidis, D. Rafferty, S. Rawlings, W. Reich, J.-P. 
de Reijer, J. Reitsma, G. A. Renting, P. Riemers, 
E. Rol, J. W. Romein, J. Roosjen, M. Rniter, 
A. Scaife, K. van der Schaaf, B. Scheers, P. Schel- 
lart, A. Schoenmakers, G. Schoonderbeek, M. Sery- 
lak, A. Shnlevski, J. Sluman, O. Smirnov, C. Sobey, 

H. Spreeuw, M. Steinmetz, C. G. M. Sterks, H.-J. 
Stiepel, K. Stnurwold, M. Tagger, Y. Tang, C. Tasse, 

I. Thomas, S. Thoudam, M. C. Toribio, B. van der 
Tol, O. Usov, M. van Veelen, A.-J. van der Veen, 

S. ter Veen, J. P. W. Verbiest, R. Vermeulen, N. Ver- 
maas, G. Vocks, C. Vogt, M. de Vos, E. van der 
Wal, R. van Weeren, H. Weggemans, P. Weltevrede, 
S. White, S. J. Wijnholds, T. Wilhelmsson, O. Wuck- 
nitz, S. Yatawatta, P. Zarka, A. Zensus, and J. van 


Zwieten (2013), LOFAR: The LOw-Freqnency AR- 
ray, Astron. Astrophys., , 556, A2, doi:10.1051/0004- 
6361/201220873. 

Vedantham, H., N. Udaya Shankar, and R. Snbrah- 
manyan (2012), Imaging the Epoch of Reioniza¬ 
tion: Limitations from Foregronnd Confusion and 
Imaging Algorithms, Astrophys. J., , 7/5, 176, doi: 
10.1088/0004-637X/745/2/176. 

Williams, C. L. (2012), Initial exploration of 21-cm cos¬ 
mology with imaging and power spectra from the 
murchison widefield array, Ph.D. thesis, Massachusetts 
Institute of Technology. 

Zaroubi, S. (2013), The epoch of reionization, in The 
First Galaxies, Astrophysics and Space Science Li¬ 
brary, vol. 396, edited by T. Wiklind, B. Mobasher, 
and V. Bromm, pp. 45-101, Springer Berlin Heidel¬ 
berg. 

Zheng, H., M. Tegmark, V. Bnza, J. S. Dillon, 
H. Gharibyan, J. Hickish, E. Kunz, A. Liu, J. Losh, 
A. Lutomirski, S. Morrison, S. Narayanan, A. Perko, 
D. Rosner, N. Sanchez, K. Schutz, S. M. Trib- 
iano, M. Valdez, H. Yang, K. Z. Adami, 1. Zelko, 
K. Zheng, R. P. Armstrong, R. F. Bradley, M. R. Dex¬ 
ter, A. Ewall-Wice, A. Magro, M. Matejek, E. Mor¬ 
gan, A. R. Neben, Q. Pan, R. F. Penna, C. M. Pe¬ 
terson, M. Su, J. Villasenor, C. L. Williams, and 
Y. Zhu (2014), Miteor: a scalable interferometer for 
precision 21 cm cosmology. Monthly Notices of the 
Royal Astronomical Society, //5(2), 1084-1103, doi: 
10.1093/mnras/stul773. 

Corresponding author: A.R. Neben, Kavli Institute for 
Astrophysics and Space Research, Massachusetts Insti¬ 
tute of Technology, Cambridge, MA 02139, USA. (abra- 
hamn@mit.edu) 

(Received_) 



NEBEN ET AL.: MEASURING MWA BEAMPATTERNS WITH ORB COMM SATELLITES 


15 



Figure 1. Simplified diagram of our beam measurement 
system and data flow. ORBCOMM satellite signals are 
received by our Reference Antenna and the AUT, each 
passing through the receiver chain described in Sec. 2 
which outputs a power spectrum with 2 kHz resolution 
between 137-138 MHz every ~200 ms. At each time step, 
a copy of predict running on our data acquisition com¬ 
puter outputs the positions and IDs of all ORBCOMM 
satellites above the horizon, while transmission frequen¬ 
cies are logged by our ORBCOMM Interface box. 
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Figure 2. Reference antenna on 2 m x 2 m ground 
screen (left) (see Section 2.3) and MWA antenna tile on 
5 m X 5 m ground screen (right) (see Section 4) de¬ 
ployed at the National Radio Astronomy Observatory- 
Green Bank. 
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Figure 3. Analysis of a typical satellite pass from 
our null experiment, in which the AUT is replaced by 
a reference-style antenna. The satellite rises out of the 
galactic background power (see Fig. 5), passes high in 
the sky, then drops below the horizon over the course of 
roughly 15 minutes (top panel). Both curves are in units 
of dB relative to the background level estimated at the 
beginning of the pass. Outside the region enclosed by ver¬ 
tical lines, the signal to background ratio is smaller than 
20 dB, meaning that the satellite signal received by each 
antenna is corrupted at the few percent and larger level 
by diffuse galactic power. The fluctuations in the ratio 
of the two antenna responses (middle panel) are typically 
at the ±0.5dB level and are due to multipath reflections 
and polarization mismatch, not receiver noise. We also 
plot the satellite zenith angle (bottom). 
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Figure 4. Results from our null experiment in which 
the AUT is another reference antenna. Beams of the EW 
(NS) oriented dipoles are shown in the left (right) col¬ 
umn. The measured AUT beampattern is plotted in dB 
relative to its boresight gain (top). These maps are in 
sine projection with North at the top and East at the 
right. Dashed circles mark 20, 40, 60, and 80 degrees 
from zenith. We also show measured and model beam- 
patterns (middle) and deviations from the model (bot¬ 
tom) on slices through the E (red) and H (blue) antenna 
planes. 
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Figure 5. Observed reference antenna background levels 
measured at the beginning or end of each satellite pass, 
plotted with a model time-dependent background com¬ 
puted from the Global Sky Model (GSM) [De Oliveira- 
Costa et al., 2008] and the model reference antenna beam- 
pattern. The data agree with the GSM model given its 
stated ±10% accuracy. 
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Figure 6. The measured MWA tile beampattern for 
the zenith pointing is plotted in dB relative to its bore- 
sight gain (top). Beams of the EW (NS) oriented an¬ 
tennas are shown in the left (right) column. The maps 
are in sine projection with North at the top and East at 
the right. Dashed circles mark 20, 40, 60, and 80 de¬ 
grees from zenith. We show measured and model beam- 
patterns (middle) and deviations from the model (bot¬ 
tom) on slices through the E (red) and H (blue) antenna 
planes. Dashed lines in the bottom panel show the model 
of Sutinjo et al. [2015] relative to the analytic model. Ver¬ 
tical dotted lines mark the model FWHM of ~23 degrees 
at 137 MHz. 
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Figure 7. Measured MWA tile power pattern for 
pointed 20 degrees West pointing. Same layout as Figure 
6, except here we plot beam slices through orthogonal 
planes through the main lobe (red and blue) which do 
not correspond to the E and H antenna planes because 
the direction of boresight is no longer in these planes. 
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Figure 8. Measured MWA tile power pattern for pointed 20 degrees East pointing. Same layout as Figure 7. 







































